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Linear theory for phonon scattering by discrete breathers in the discrete nonlinear 
Schrodinger equation using the transfer matrix approach is presented. Transmission and 
reflection coefficients are obtained as a function of the wave vector of the input phonon. 
The occurrence of a nonzero transmission, which in fact becomes perfect for a symmetric 
breather, is shown to be connected with localized eigenmodes thresholds. In the weak- 
', coupling limit, perfect reflection are shown to exist, which requires two scattering channels. 

■ A necessary condition for a system to have a perfect reflection is also considered in a general 
' context. 

^ ■ I. INTRODUCTION 

(N 

, Breathers, -which also called oscillating solitons, are time periodic and spatially localized excitations ubiq- 

^ , uitously found on translationally invariant coupled nonlinear lattice 0-0]. It -was introduced first in the 

■ context of the sine-Gordon equation in 1974 It requires practically no activation energy and thus bridge 
I the gap bet-ween the highly nonlinear modes and the linear phonon modes. Such long-lived spatially localized 

excitation has fundamental significance in the dynamics of spatially extended systems. It has been found 
that the discreteness of a system plays an essential role on the stability of breathers in several systems ||] . 
. It introduces a cut-off frequency to the system so that the resonance bef-ween the phonon band and the 

' frequency of the discrete breather can be avoided. 

Ho-wever, breathers gro-w and decay due to their interaction -with the environment. Some of such interac- 
tions -which has been studied include the electron-breather interaction Q and the interaction of breathers 
-with extended defect 0. The discrete breather can play a role of the scattering center affecting energy 
transport by scattering, absorbing or radiating phonons. Since the phonon is a major form of radiations for 
Q ^ small amplitude excitations in the nonlinear lattice, it is important to characterize such breather-phonon 

^ . interaction 

Recently it has been found that scattering properties of the breather are related to the structure of the 
- ^ , internal modes of the breather itself j|] . They studied the phonon scattering by a kink in the localized 

' equilibria of the nearest neighbor interacting chain and showed that the scattering properties undergo a 

^ I drastic change, involving a perfect transmission, as the localized eigenmode threshold is passed. It -was 

conjectured that this result on the static localized structure could be extended to a -wide class of systems 
-with time-dependent localized structures such as breathers. It has also been confirmed in the case of time 
periodic discrete breather in the Hamiltonian nearest neighbor chain both numerically and using the extended 
Levinson theorem [^j. 

The purpose of this paper is to present the analysis of phonon scattering by breathers in the discrete 
nonlinear Schrodinger (DNLS) equation using the transfer matrix method. We proved analytically that at 
the localized eigenmode threshold the transmission at the phonon band edge, such that the wave vector 
fc = 0, is nonzero, -which in fact becomes perfect in the case of symmetric breather. This result is also 
confirmed by numerical calculations -with transfer matrices. The crucial difference of the DNLS equation 
from the localized equilibria of the nearest neighbour chain is that it allo-ws f-wo scattering channels H . This 
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allows occurrence of the perfect reflection whose existence is shown in detail in the weak-coupling limit. A 
necessary condition for a system to have a perfect reflection is obtained in comparison with the case of the 
local equilibria in the classical Hamiltonian systems. We end with discussions on extensions of our results 
to more general systems including the time-periodic discrete breathers in the Hamiltonian chain. 



II. THE TRANSFER MATRIX FORMULATION 



There are many ways to discretize the nonlinear Schrodinger equation. One of the widely studied version, 
called the DNLS equation, reads, 

= i{\''Pn\'^1p7i + eii^n-l + i'n+l)}, (1) 

in which ipn is the amplitude at site n and e is the coupling strength. If one consider a solution of the form 
V'n = One*"* for Eq.(0), a„'s satisfy 

wa„ = -(- e(a„_i + a„+i). (2) 

Without loss of generality, we could choose w > 0. We further consider only the 1-site breather, i.e. only one 
site is excited on the lattice. Hence for zero-coupling the breather solution centered at the origin is given by 
ao = y/uj and a„ = for n 0. 

In the weak-coupling limit, i.e. if e/w is small, we obtain from Eq.(|2|) oq = v^(l ~ (^Z"^)^ + ' ' ^±1 = 
■^/Lj{e/uj +•■■), ■ ■ ■ etc. This provides an exponentially localized amplitude profile of the discrete breather 
which can be continued from a zero-coupling solution pO| ]. The reflection symmetry of the breather with 
respect to its center will be shown to be important in connection with perfect transmission. 

The linearization of Eq.(P around the breather solution with 4>n = 4'n — ane^^^ yields, 

(/)„ = i{2al(bn + ay^^'r„ + e{(bn+i + (3) 

Because of the complex conjugate operation on 0„, at least two frequency components are necessary for the 
solution of Eq.(^). Therefore, the DNLS equation can allow two scattering channels different from the case 
of localized equilibria with a single scattering channel If we assume </>„ = &„e*'"+^-'* -I- c„e*^"~^^*, with 
bn and c„ complex, we obtained the following equations: 

' + n-2al \ f bn\ ( al\ ( hl\ _ fb„-i+b„+ 



With = {bn, c„, 6* , c* ), the above equation can be written in a recursive form as 

^" ^=mJ\+'), (5) 



in which M„ is a function of a„, e, lo and Q. The matrix Af„ is called a transfer matrix, which maps a pair 
of linearized coordinates on one site to another site Applying Eq.(|^) successively from the site to 

the site —N, one has 

^=m(IZ^A, (6) 



with M = M_Ar • M_Ar+i • • • Mn^i ■ Mm- 

Far from the center of the breather, a„ decays to zero, so that Eq.(|^) becomes decoupled, with solutions 
of the form &„ = ^^kin _ ^ik^n ^ jj-^ which cosfcfc = (tj + il)/2e and cos fee ~ {lo — r2)/2e. The solutions 

of bn and c„ define two possible scattering channels. Each channel can be either a traveling (phonon) mode 
with a real wave vector k or an exponentially growing (decaying) mode with an imaginary k depending upon 
the values of \u) ± Q\/2e. Since e and O occurs only in the form of ratios e/uj and fi/o;, we can set lo — 1 and 
regard e and 51 as dimensionless variables. The boundaries of the phonon bands are given by w = il ± 2e 
and — r2 ± 2e. These phonon bands divided the parameter space of e and f2 into three regions according to 



the number of traveling channels as shown in Fig.Vl 
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III. THE SCATTERING PROBLEM 



In a Levison theorem, the bound state property of a system is related to its scattering property p^ . In 
what foUows, the property of scattering in region-1 was related to the property of localized mode in region-0. 

Firstly, we consider the scattering problem in region-1. In the standard scattering setup, a phonon is 
incident with normalized amplitude from the left side of the breather located at the origin. We assume 

> 0. The case with 51 < can be handled similarly. 

The question asked in a scattering problem is its transmission and reflection coefficients. It can be shown 
from the DNLS equation that 



(7) 



in which < A > is the time average of A. This flux-conservation relation makes the one-channel scattering 
process elastic, i.e., 



(8) 



in which t and r are transmission and reflection coefficients respectively. 

Suppose the channel c„ is a traveling mode with kc real, while kb — iKbinb > 0) is imaginary. The 
asymptotic solutions are of the form 



bn 



= b+t 



if; 



if n 



(9) 



in which the amplitude of the incident phonon on the channel c„ is set to unity. The amplitudes for the 
transmitted and reflected phonons, and , fully characterize the one-channel phonon scattering by a 
symmetric breather. Since Eq.(^ is not covariant under multiplication of a arbitrary phase factor e'^ on 5„ 
and c„, it may appear that the scattering properties depend on the phase of the input phonon. However, if 
bn and c„ are solutions of Eq.(|4|), so are ibn and —icn- Therefore, transmission and reflection are independent 
of the phase of the input phonon in the one-channel scattering case. 

Inserting Eq.(||) into Eq.(||), we obtain 8 linear equations with 8 variables to solve for, i.e. 



-ikaN 



„ikaN 



-ik^N 



5-e-«5(w+i) 



ka{N+l) 

I ^ 

y ik^(N+l) _|_ g-*g-ifce(W+l) j 



= M 



g+gjA;,(W+l) 
5+'*g-Kb(A'+l) 
g+*g-ifc<:(Ar-|-l) 



bt* 



-KtN 



l^+*g—ikaN 



(10) 



For kc — 0, 1 + can be treated as one complex variable. The resulting 8 homogeneous equations have a 
non-trivial solution only on their null space. 

On the other hand, if both channels are non-traveling modes with both kb — iKb^Kb > 0) and k^ = ittci^c > 
0) imaginary. After discarding all exponentially growing parts, the asymptotic solutions read. 



bi_e 
Ci e 
b+t 



if n 
if n 



(11) 



oo. 



The localized eigenmode is given by a solution of 
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V 



(12) 



Non-trivial solutions for this set of equations exist only on its null space. At fee = the null spaces of Eq.(]lO|) 
and Eq.(jl^) coincide. Therefore, at the band edge of the channel c„, fee — 0, where a localized mode start 
to appear, we obtain non-zero phonon transmission, c+. 

We can further demonstrate that for symmetric breathers non-zero transmission at fcc = at the threshold 
for localized mode generation becomes perfect one: Consider a phonon of unity amplitude incident from the 
left side of the breather at a small positive wave vector fee- Near = 0, this scattering configuration can be 
extended continuously to a small negative value of since the governing equations are smooth functions of 
kf.. Negating /cc changes the direction of traveling phonons. The configuration becomes an outgoing phonon 
of unity amplitude and an incoming phonon of amplitude r{—kc) on the left side of the breather, and an 
incoming phonon with amplitude t{—kc) on the right side as shown in Fig.^. For a symmetric breather, 
these incoming phonons with a wave vector —kc from the right side of the breather generate outgoing phonons 
of amplitude fcc) + P{—kc) and 2r{—kc)t{—kc) from the left and right sides of the breather, respectively. 
The situation in Fig.||b and Fig.Hc must coincide for kc = 0. Therefore, we obtain r^(0) -I- t^(0) = 1 and 
2r{Q)t{0) = 0, and conclude that there must be either perfect transmission or perfect reflection at kc = 0. 
This is a special property of the symmetric breather since we assumed that the transmission and reflection 
coefficients for a phonon coming from the left side of the breather are the same as those from the right side 
of the breather. 

In general, a scaterring process can be formulated as a set of homogeneous linear equations for the Fourier 
components of the state of the system at each site and casted into a transfer matrix form if the couplings 
between sites are local and linear. Let K„ be the state vector at site n and Af„ be the transfer matrix 
connecting {Y„+i,Yn)'^ and (F„,y„_i)^ as in Eq.(||). If the dimension of y„ is one, the whole vector F„ 
represents the traveling phonon in the scattering process. If the incoming phonons from the left side of the 
breather are perfectly reflected, the Fjv should be identical to zero as — + oo. Because of the homogeneity 
of Eq.(^, y_Ar should also be zero, which implies no non-trivial configuration of perfect reflection is possible. 
If the dimension of y„ is greater than one and some components represent the exponentially decaying parts, 
Yn need not to be zero in the N oo limit in spite of the perfect reflection. It becomes possible to 
match non-zero K-at, which represents the input and reflection of phonons, with Y/v by the transfer matrix. 
Therefore, the existence of additional channels with exponentially decaying parts is essential for a perfect 
reflection to occur in a linear scattering problem. Hence, there is no perfect reflection in two-channel phonon 
scattering of the DNLS equation, i.e. in region-2 of Fig. VI, since all components of Yn represent traveling 
phonons. This provides a necessary condition for a linear-system to allow a perfect reflection. 

Our derivation does not depend on the specific form of the transfer matrix AI. The result is entirely 
based on the fact that asymptotic solutions which are composed of traveling phonon parts and exponentially 
growing or decaying parts far from the center of the breather can be connected to each other by a transfer 
matrix. Therefore, the same procedure can be applied to the breather-phonon scattering problem in the 
nearest neighbor Hamiltonian chain to obtain a similar proof [||. A major difference is that the number of 
channels in the latter case should be extended to infinity for each frequency component f2 -I- ojf,, in which f2 
is the phonon frequency and Wf, is the breather frequency, in the case of Hamiltonian breather. 



IV. PERTURBATIVE AND NUMERICAL CALCULATIONS 



The occurrence of the perfect reflection in the weak-coupling limit can be shown by a perturbation method. 
Retaining only terms of order e/ui, we obtain the transfer matrix M « M^i-Mq-Mi; other M„'s are discarded 
because they affect only the phase factors in the order e/co. The transmission coefficient thus obtained is 
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(13) 



The breather becomes opaque to the fourth power of the couphng strength. A perfect reflection occurs at 
k(. = 7r/2 in the weak-couphng hmit. Numerical calculation, which takes into account of several higher order 
terms of e/^jJ, agrees with the perturbative result as shown in Fig.^. Perfect phonon reflection by breathers 
in the DNLS equation is a property which persists for weak-coupling, in contrast to what happens in the 
Hamiltonian nearest neighbor chain where the perfect reflection exist only for a positive interval of coupling 
strength e for a given Wf,. 

In the case of phonon scattering by a local equilibria kink in the nearest neighbor Hamiltonian chain, there 
is no perfect reflection However, there occurs a perfect reflection in the Hamiltonian discrete breather 

. The critical difference between the case of a local equilibria kink and the case of Hamiltonian or the 
DNLS breathers is the minimum number of frequencies needed for describing linearized motion at each site. 
In the case of the local equilibria, one Fourier component measuring the deviation from the equilibrium at 
each site is sufficient to describe the motion of each site because the kink in the local equilibria is a time 
independent defect. On the other hand, the discrete breather in the Hamiltonian system has an intrinsic 
time scale Tf, which is the period of the breather. The intrinsic frequency — ^ix jT-h generates an infinite 
number of frequencies + mcjb (jn is an integer ) mixing with the phonon frequency by nonlinear terms 
in the equation. For the DNLS case, only one more frequency is generated due to the complex conjugate 
operation in Eq.(^). 

With the transfer matrix method, the transmission and reflection coefficients can be calculated numerically 
in the same way as in the case of the local equilibria The transfer matrix, M, is obtained numerically 
and the scattering coefficients are solved from Eq.([l^). The number of sites in the scattering setup, N , in 
Eq.(|lO|) is chosen to allow a decay of the amplitude of the breather within the error bound, 10^^, in our 
calculation. For 0.01 < e < 0.4, iV needs not to exceed 5 to maintain this bound. 

The phase diagram for the localized eigenmodes and the transmission for various values of e are obtained 
numerically in Fig.^ and Fig.^. For small coupling, there is no perfect transmission. A perfect transmission 
is generated from /cc = at e = 0.262 where a symmetric localized eigenmode fl^ is generated from the 
lower edge of the phonon band. The value of kc where the perfect transmission occurs increases from zero 
until it reaches a maximum value. Then the maximal wave vector decreases and becomes zero again where 
another localized eigenmode, which is anti-symmetric at this time | p3| , is generated. These numerical results 
are consistent with predictions of the linear scattering analysis. 

The return of fee, where the perfect transmission occurs, to zero after its initial creation and deviation 
from zero is due to the multi-site spreading of the breather envelope. For a one-site breather, if only Mq is 
retained in calculating M, the location of the perfect transmission does not return to /cc = again; which 
is consistent with the fact that there is no generation of the anti-symmetric localized eigenmode. In other 
words, the generation of a perfect transmission and a symmetric localized eigenmode is due to the single- 
site defect property of the discrete breather, whereas the annihilation of the perfect transmission and the 
generation of an anti-symmetric localized eigenmode is due to the multi-site structure of the breather. 

The perfect reflection at nonzero k^ is another special feature of the one-channel phonon scattering case. 
It is important because the breather can be regarded as a perfect reflecting barrier for the phonon, which 
opens the possibility of trapping phonons between two breathers. Numerical results show that a perfect 
reflection occurs at a special value of k^ if the coupling is not too strong. The value of kc where the perfect 
reflection occurs changes from 7r/2 to as e increases. If e reaches a critical value of 0.36, the perfect reflection 
disappears through fee = as depicted in FigJ|. 

V. CONCLUSION 

In conclusion, phonon scattering by discrete breathers in the DNLS equation is formulated by the transfer 
matrix method and transmissions were obtained in the one-channel scattering case. The nonzero transmission 
of phonons through the discrete breather at /c = in the one-channel phonon scattering occurs at the localized 
eigenmode threshold is proven analytically to be perfect for the symmetric discrete breather. Numerically, 
the perfect transmission was shown to occur at /c = whenever the symmetric and antisymmetric localized 
eigenmodes were generated from the lower edge of the phonon band. We obtained explicitly the formula for 
the phonon transmission in the lowest order of the coupling strength. A perfect reflection was shown to exist 
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at fc = 7r/2 in the small coupling limit. The location of the perfect reflection was traced numerically until it 
disappeared through fee = at a critical value of the coupling strength. 

The solution matching method using the transfer matrices has been proved to be an efficient tool for 

investigating scattering properties. Furthermore, the DNLS system provided a good intermediate example 
with two scattering channels bridging the gap between the case of a kink in the local equilibria with a single 
scattering channel and the discrete breather in the classical Hamiltonian system with an infinity of scattering 
channels. The methods used for the breather in the DNLS system can be extended to study the one-channel 
scattering of the nearest neighbor Hamiltonian chain. The multi-channel scattering problem is important 
in the context of stability, e.g., a two-channel scattering problem where the phonons of the second channel 
are incident from one side of the breather. Detail study of phonon-induced instability will be the subject of 
on-going investigations. 
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FIG. 1. The eigenspectrum of linearized modes of the DNLS equation as a function of e. The phonon bands are 
shaded. The number in each region of the diagram denotes the number of the traveUng channels. 
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FIG. 2. Continuation of the one-channel phonon-breather scattering configuration from a small positive value of 
kc to a small negative one. 
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FIG. 3. Phonon transmission through the breather in the DNLS equation in the weak-couphng limit. The sohd 
hne is the perturbative result, while diamonds denote the result of numerical calculations using the transfer matrix 
method. 
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FIG. 4. The blowup of Fig.^ with two branches of locahzed eigenmodes coming out of the phonon band. 
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FIG. 5. One-channel phonon transmissions through discrete breathers in the DNLS equation for w = 1. (a)-(b) 
Near e = 0.262 a perfect transmission is generated, (c)-(d) Near e = 0.315, the perfect transmission is annihilated. 
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FIG. 6. One-channel phonon transmissions through breathers in the DNLS equation for w = 1. Near e = 0.36, (a) 
the perfect reflection occurs at or tt; (b) as the couphng strength is increased it is annihilated at A; = 0. (c)The 
wave vector for the perfect reflection as a function of the coupling strength e. 
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